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Abstract 



With its three independent parameters, the ABCD matrix serves 
' (— I '■ as the beam transfer matrix in optics. If it is transformed to an 

Qh' equi-diagonal form, the matrix has only two independent parameters 

fH \ determined by optical devices. It is shown that this two-parameter 

mathematical device contains enough information to describe the ba- 
sic space-time symmetry of elementary particles. If its trace is smaller 
than two, this matrix can represent the internal space-time symmetry 
of massive particles. If equal to two, the matrix is of massless parti- 
^ ' cles. If the trace is greater than two, this matrix describes imaginary- 

mass particles. This matrix speaks Einstein's language for space-time 
structure of elementary particles. As for the optical devices, the laser 
ff^ ■ cavity and the multilayer system are discussed as illustrative physical 

p^i ■ examples. 

1 Introduction 

^ • The ABCD matrix is a two-by-two matrix with real elements, and its de- 

. terminant is one. There are therefore three independent parameters. These 

elements are determined by optical materials and how they are arranged. 
The purpose of this note is to explore its mathematical properties which can 
address more fundamental issues in physics. 

First of all, the trace of this matrix could be less than two, equal to two, 
or greater than two. We are interested in what physical conclusions we can 
derive from these numbers. 
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In order to bring the ABCD matrix to the form which will address some 
fundamental issues in physics, we should first transform it into the equi- 
diagonal form where the two-diagonal elements are equal to each other [11 [2]. 
We can achieve this goal by a similarity transformation with a one-parameter 
matrix. This transformation does not change the trace, and the resulting 
equi-diagonal matrix has two independent parameters. 

We shall call this equi-diagonal matrix the core of the ABCD matrix, 
and use the notation [ABCD]. This matrix cannot always be diagonalized. 
This creates a non-trivial problem. We shall examine how optical devices, 
especially periodic systems, can lead us to a better understanding of the 
problem. For this purpose, we discuss laser cavities and multilayer systems 
in detail. 

If the trace is less than two, the core can be written as 



The diagonal elements are equal and smaller than one. 

If the trace is greater than two, the [ABCD] matrix takes the form 



Here again the diagonal elements are equal, but they are greater than one. 
If the trace is equal to two, the core matrix becomes 



This matrix also has the same diagonal element, and they are equal to one. 

This triangular matrix of Eq.Q cannot be diagonalized. The core ma- 
trices of Eq.([T]) and Eq.([2]) can be diagonalized, but not by rotation alone. 
These mathematical subtleties are not well known. The purpose of this re- 
port is to show how much physics we can extract from these mathematical 
details. 

The mathematics of group theory allows us to write down a four-by-four 
Lorentz-transformation matrix for every two-by-two matrix discussed in this 
paper. In this way, the three matrices given in Eq.([T]), and Eq.([2]), and 
Eq.Q lead to the internal space-time symmetries of elementary particles. 
They respectively correspond to the symmetries of massive, imaginary-mass, 
and massless particles respectively [3]. 




(1) 






(3) 
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In Sec. m we write the [ABCD] matrix in terms of the its generators, and 
decompose it to three matrices in the form of similarity transformation. It 
is noted that there is another mathematical device known as the Bargmann 
decomposition [1]. This decomposition is not a similarity transformation, 
but it can play other important roles in understanding the ABCD matrix. 

In Sec. El we discuss a laser cavity as a physical example of the math- 
ematical details of the [ABCD] matrix. In Sec. lU a multilayer system is 
discussed as a physical example leading to the desired similarity transforma- 
tion. The Bargmann decomposition plays the key role in this problem. In 
Sec. [5l we shall discuss the space-time symmetries implied by the properties 
of the [ABCD] matrix. The Bargmann decomposition, as well the similarity 
decompositions, is explained in terms of the Lorentz transformations which 
leave the momentum of a given particle invariant. Thus, these transforma- 
tions are applicable to internal space-time symmetries. 

2 Decomposition of the ABCD Matrix 

We are interested in writing the three different forms of the core matrix in 
one expression. 



where the parameters x and y are determined by the optical materials and 
how they are arranged. The exponent of this matrix becomes 






x-y 




(5) 



If X > ?/, the matrix becomes 




exp{r]/2) 





which leads to the core matrix of Eq.([T]) with 





(7) 
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The core matrix [ABCD] can be written as a similarity transformation 

[ABCD] = B{r])R{e)B{-7]) (8) 

with 

Biv) = [ 

/cos(0/2) -sin(e/2)\ 
^ ' \sm{e/2) cos(^/2) r 

where 7 is now replaced by the rotation angle 6. R{9) is a rotation matrix, 
and -8(77) is a squeeze matrix. 
If X < the exponent becomes 



If exp (r//2) \ 

2Vexp(-r//2) j 

leading to the core matrix of Eq. ([2]) , with 



(10) 



7 = vi/^ - ^^ 



e' = (11) 

The [y4i?CD] matrix can now be decomposed into a similarity transformation 

[ABCD] = B{r^)S{-X)B{~r^), (12) 

with 

/cosh(A/2) sinh(A/2)\ 
•^^^^ ~ Uin(A/2) cosh(A/2)j' ^^'^^ 

where 7 is replaced by the boost parameter A. The matrix B{rj) takes the 
diagonal form given in Eq.(|8]) with 7] defined in Eq. ffTT]) . S'(A) is a squeeze 
matrix. 

U X = y, the exponent becomes 



1 —X 

1 



(14) 



with X = y = J. 
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We now have combined three different expressions for the core of the 
ABCD matrix into one exponential form of Eq.(jl]). This form can be de- 
composed into three matrices constituting a similarity transformation. We 
shall call this the 'Wigner decomposition" for the reasons given in Sec. |5l 

There is another form of decomposition known as the Bargmann decom- 
position [1], which states that the core of the ABCD matrix can be written 
as 

[ABCD] = R{a)S{-2x)R{a), (15) 

where the forms of the rotation matrix R and the squeeze matrix S are given 
in Eq.([8]) and Eq.(|9]) respectively. If we carry out the matrix multiplication, 
the [ABCD] matrix becomes 



(cosh x) cos a — sinh x — (cosh x) sin a 

— sinh X + (cosh x) sin a (cosh x) cos a 




This matrix also has two independent parameters a and x- We can write 
these parameters in terms of 7 and rj by comparing the matrix elements. For 
instance, if x > the diagonal elements lead to 

cos(^/2) = (coshx) cosa. (17) 

The off-diagonal elements lead to 

^ {coshx) sin a + sinh x 
(coshx) sin a — sinh x 

As for physical applications, let us consider periodic systems, such as laser 
cavities and multilayer systems. The exponential form given in Eq.(|l]) tells 
us that it is a matter of replacing the 7 parameter by A^7 for N repeated 
applications [2]. Let us see some examples. 



3 Laser Cavities 

As the first example, let us consider the laser cavity consisting of two identical 
concave mirrors separated by a distance d. Then the ABCD matrix for a 
round trip of one beam is 



1 0\fl d\f 1 0\ fl d 
-2/R 1) [0 1) \ -2/R 1) [0 1 



(19) 
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where the matrices 

are the mirror and translation matrices respectively. The parameters R and 
d are the radius of the mirror and the mirror separation respectively. This 
form is quite familiar to us from the laser literature O El 13 |8] . 

However, the crucial question is what happens when this process is re- 
peated. We are thus led to the question of whether the chain of matrices in 
Eq. (fTO!) can be brought to an equi-diagonal form and eventually to a form of 
the Wigner decomposition. We are interested in finding the core of Eq. (fT9|l . 
For his purpose, we rewrite the matrix of Eq. (fT9|) as 

1 -rf/2\/l d/2\ f 1 OWl d/2\^ 



1 j VO 1 ) \-2/R IJ \0 1 

In this way, we translate the system by —d/2 using a translation matrix 
pven in Eq.f l20p . and write the ABCD matrix of Eq. ffTOj) as 



1 -d/2\ 



1-d/R d-d^/2R\ 



1 d/2 
1 



,0 I ) [\ -2/R 1-d/R ) 
We are thus led to concentrate on the matrix in the middle 



(22) 



1-d/R d-dy2R\ , , 

-2/R 1-d/R )' 



which can be written as 



y/d \fl-d/R 1 - d/2R\ fl/Vd 
l/Vd)\-2d/R 1-d/R J [ Vd 



(24) 



It is then possible to decompose the ABCD matrix into 

E E-\ (25) 



with 

C 



1-d/R l-d/2R\ 
-2d/R 1 - d/R ) 
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The C matrix now contains only dimensionless numbers, and it can be written 

as 

/ cos(7/2) e''sin(7/2)\ 

l-e-^sin(7/2) cos(7/2) J' ^ 



with 



cos(7/2) = 1 - |, 



Here both d and R are positive, and the restriction on them is that d be 
smaller than 2R. This is the stability condition frequently mentioned in the 
literature [H [7j . 

Thus, the [ABCD] core matrix is C^, and takes the form 

\A^>rn^-^ e''sin(7)\ , . 

Ue-''sin(7) cos(7) j' ^^^^ 

and the similarity transformation which connects this core matrix with the 
original ABCD matrix of Eq.( |T9i) is the matrix E given in Eq.( l26l) . 



4 Multilayer Optics 

We consider an optical beam going through a periodic medium with two 
different refractive indexes. If the beam traveling in the first medium hits 
the second medium, it is partially transmitted and partially reflected. In 
order to maintain the continuity in the pointing picture, we normalize the 
electric field as 

iJp'' = exp (±ikiz — ut) 



E^^^ = exp {±.ik2Z - cot) (30) 



^2 

for the optical beams in the first and second media respectively. The super- 
script (+) and (— ) are for the incoming and reflected rays respectively. 
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These two optical rays are related by the two- by-two ABCD matrix, 
according to 

A B 
C D, 





(31) 



Of course the elements of this matrix are determined by transmission coeffi- 
cients as well as the phase shifts the beams experience while going through 
the media [9]. 

When the beam goes through the first medium to the second, we may 
use the boundary matrix[ll] 

n(rr\ - ( cosh(a/2) sinh(a/2) \ , . 

'^^'^^ ~ Uinh(a/2) cosh (a/2) J ' ^"^^^ 

where the parameter a is determined by the refiection and the transmission 
coefficients [S], HI [TTj. Then the boundary matrix for the beam going from 
the second medium should be Q{—a). 

In addition, we have to consider the phase shifts the beams have to go 
through. When the beam goes trough the first media, we can use the phase- 
shift matrix 

^(^i)=( e^'/O' ^^^^ 
and a similar expression for P {82) for the second medium. 

We are thus led to consider one complete cycle starting from the midpoint 
of the second medium, and write 

P (<52/2) Q{a)P (5i) Q{-a)P {82/2) . (34) 

If multiplied into one matrix, is this matrix equi-diagonal to accept the 
Wigner and Bargmann decompositions? Another question is whether the 
matrices in the above expression can be converted into matrices with real 
element. 

In order to answer the second question, let us consider the similarity 
transformation 

Ci P{5)Q{a) (35) 

with 

--;^(M)- (-) 

This transformation leads to 

R{5)Q{a), (37) 
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where 

•cos(5/2) -sm(5/2)\ 



^^^^ - [ sin(5/2) cos(5/2) j ' ^^^^ 

This notation is consistent with the rotation matrices used in Sec. El 
Let us make another similarity transformation with 

This changes Q{(j) into B{a) without changing R{6), where 

5(-)=('o' (40) 

again consistent with the B{ri) matrix used in Sec. [2j 
Thus the net similarity transformation matrix is (TU] 

C = C2C1 = _e-^^/4 e"^''/^ j ' '^'^^^ 

with 



^ ^ ~ ~ ( in/4, ) • (42) 



1 / e~^^^^ _g«'r/4 

If we apply this similarity transformation to the long matrix chain of 
Eq.(l34l). it becomes another chain 



M = R (52/2) B{a)R (5i) B{-a)R {62/2) , (43) 

where all the matrices are real. 

Let us now address the main question of whether this matrix chain can 
be brought to one equi-diagonal matrix. We note first that the three middle 
matrices can be written in a familiar form 

M = B{a)R{6i)B{-cj) 

cos((5i/2) -e'^sin(5i/2)\ 
e-"sin(5i/2) cos(5i/2) J' 

However, due to the rotation matrix R{62/2) at the beginning and at the 
end of Eq. fH3l) . it is not clear whether the entire chain can be written as a 
similarity transformation. 
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In order to resolve this issue, let us write Eg. fl44p as a Bargmann decom- 
position 

Ria)Si-2x)Ria), (45) 

with its explicit expression given in Eq. (fT6l) . The parameters a and x ^''^^ 
related to a and by 

cos(5i/2) = (coshx) cosa, 

(46) 



2^ (coshx) sin a + sinh x 



{coshx) sin a — sinh x 

It is now clear that the entire chain of Eq.f lS^ can be written as another 
Bargmann decomposition 

M = R{a + 62/2)S{-2x)R{a + 62/2). (47) 

Finally, if the trace of this matrix is less than two, it can be converted to a 
Wigner decomposition 

M = B{r])R{9)B{-7]) (48) 

with 

cos(^^/2) = (coshx) cos(a + 52/2), 

h -v 

(49) 



^2r, _ (cosh x) sm{a + 62/2) + sinh x 



(coshx) sin(a + 62/2) — sinhx 

Similar expressions can be derived for other values of the trace [12] . 

It is interesting to note that the Bargmann decomposition plays an essen- 
tial role in bringing the chain of five matrices given in Eq. fH^ to the Wigner 
decomposition of Eq. fHS]) consisting three matrices. We shall see what other 
function the Bargmann decomposition can perform in Sec. [51 



5 Space-time Symmetries 

These properties are applicable to many other branches of physics. For in- 
stance, one of the persisting problems is the internal space-time symmetry 
of elementary particles in Einstein's Lorentz-covariant world. 
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The mathematics of group theory allows us to translate the rotation and 
squeeze matrices of Eq. and Eq. into the following four-by-four matrices 
respectively. 



R{9) 



S{X) 



^ COS 9 


n 
u 


sin^ 0\ 





1 







— sin 6' 





cos^ 


V 







ij 


/ cosh A 








sinh A \ 





1 














1 





\ sinh A 








cosh A / 


/I 









1 













cosh 


V 


sinh?7 


Vo 


sinh 


V 


cosh T] J 



(50) 



They are applicable to the Minkowskian four- vector (x, y, z, t). The R{9) 
matrix performs a rotation around the y axis, and >S'(A) is for Lorentz boosts 
along the x axis. The B{ri) matrix boosts the system along the z direction. 

Together with a rotation matrix around z axis 



Z(0) 



/ COS( 

sin( 


V 



— sm I 

coscj) 





0\ 



1 

ly 



(51) 



they constitute Wigner's little groups dictating internal space-time symme- 
tries of massive and imaginary- mass particles [3]. The triangular matrix of 
Eq.(|3]) leads to the little group for massless particles. The little groups are 
the subgroups of the Lorentz group whose transformations leave the four- 
momentum of a given particle invariant. 

Let us go back to Eq.([T]) which, according to Eq.®, can be decomposed 
to a similarity transformation 



W{ri,9) = B{ri)R{9)B{-ri). 



(52) 



We can write this decomposition with the four-by four matrices given in in 
Eq.(ISnD. 
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Let us then consider a massive particle moving along the z direction with 
the velocity parameter v/c = tanh?7, and its four-momentum 

(0, 0, m sinhi], m cosh r/), (53) 

where m is the mass of the particle. 

We can boost this particle using the boost matrix B^—t]), which is the 
inverse of the four- by-four matrix given in Eg. (1501) . The particle becomes at 
rest, with its four- momentum 

(0,0, 0,m), (54) 

and with zero velocity. The rotation matrix R{6) rotates this particle with- 
out changing its momentum. During this process, the particle changes the 
direction of its spin. Finally, B{ri) boosts the particle and restores its mo- 
mentum, as is illustrated in Fig. [H Thus, it is appropriate to call the form of 
Eq. fl52l) the Wigner decomposition. In this way, the four-by-four expression 
for Eq.([H]) changes the internal space-time structure of the particle. 

The essential function of the Wigner decomposition is to provide the 
subgroup of the Lorentz group which will leave the given four-momentum of 
a particle invariant [3]. 

The Bargmann decomposition also provides momentum-preserving trans- 
formations. The decomposition consists of a rotation, a boost, and another 
rotation as illustrated in Fig. [TJ It is interesting to note that the Wigner de- 
composition and the Bargmann decomposition can serve the same purpose 
of providing a Lorentz transformation which leaves the momentum invariant. 

One key question from this table is what happens to the 0(3)-like little 
group when the particle momentum becomes infinity or its mass becomes 
zero. The question is whether the little group for a massive particles becomes 
that for a massless particle. The answer to this question is Yes, but this 
issue had a stormy history before this definitive answer [T3]. Indeed, when rj 
becomes infinity, the four-by-four form of Eq (l52l) becomes 



/I 





-7 


7 \ 





1 








7 





1 - 772 


772 
1+772; 


V7 





-772 



(55) 



When applied to the momentum of a massless particle moving in the negative 
z direction with 

(0,0,p,p), (56) 
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Figure 1: Wigner decomposition (left) and Bargmann decomposition (right). 
These figures illustrate momentum preserving transformations. In the 
Wigner transformation, a massive particle is brought to its rest frame. It 
can be rotated while the momentum remains the same. This particle is then 
boosted back to the frame with its original momentum. In the Bargmann 
decomposition, the momentum is rotated, boosted, and rotated to its original 
position. 

This matrix leaves the above four-momentum invariant, but it performs a 
gauge transformation [I3]. This aspect of Wigner's little group is illustrated 
in Table [H 
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Table 1: Massive and massless particles in one package. Einstein unified the 
energy-momentum relation for slow (massive) and fast (massless) particles 
with one Lorentz-covariant formula. Likewise, Wigner's little group unifies 
the internal space-time symmetries of particles. When boosted, the spin 
along the direction remains invariant. This is called the Lorentz- invariant 
helicity. The spins along the transverse directions collapse into a gauge degree 
of freedom in the infinite- momentum or zero- mass limit [131 [E] • 





Massive 


Lorentz 


Massless 




Slow 


Covariance 


Fast 


Energy- 


E = 


Einstein's 




Momentum 


/2m 


E= [p2 + ^2]l/2 


E = p 


Internal 


Ss 


Wigner's 


S3 


Symmetries 


Si, 5*2 


Little Group 


Gauge Trans. 



Concluding Remarks 

In this report, we have discussed some properties of the ABCD matrix which 
serves as the standard research tool in ray optics. This two-by-two matrix 
four real elements, but only three independent parameters if the determinant 
of the matrix is constrained to be one. 

If the determinant is restricted to be one, the ABCD matrix has three 
independent parameters. It was noted that this matrix can be written as 
a similarity transformation of a core matrix with two independent parame- 
ters. Two physical examples are given to illustrate how these parameters are 
determined from optical devices. 

It is remarkable that this two-parameter matrix contains enough infor- 
mation to describe the internal space-time structure of elementary particles. 
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